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ABSTRACT 

We combine data from The H I Nearby Galaxy Survey and the GALEX Nearby Galaxy Survey to 
study the relationship between atomic hydrogen (H I) and far-ultraviolet (FUV) emission outside the 
optical radius (r25) in 17 spiral and 5 dwarf galaxies. In this regime, H I is likely to represent most 
of the ISM and FUV emission to trace recent star formation with little bias due to extinction, so 
that the two quantities closely trace the underlying relationship between gas and star formation rate 
(SFR). The azimuthally averaged H I and FUV intensities both decline with increasing radius in this 
regime, with the scale length of the FUV profile typically half that of the H I profile. Despite the 
mismatch in profiles, there is a significant spatial correlation (at 15" resolution) between local FUV 
and H I intensities; near r25 this correlation is quite strong, in fact stronger than anywhere inside r25 
(where H I is not a good tracer for the bulk of the ISM), and shows a decline towards larger radii. 
The star formation efficiency (SEE) — defined as the ratio of FUV/H I and thus the inverse of the gas 
depletion time — decreases with galactocentric radius across the outer disks, though much shallower 
than across the optical disks. On average, we find the gas depletion times to be well above a Hubble 
time (~ lO^^yr). We observe a clear relationship between FUV/H I and HI column in the outer 
disks, with the SEE increasing with increasing H I column. Despite observing systematic variations 
in FUV/H I, we find no clear evidence for step-function type star formation thresholds, though we 
emphasize that it may not be realistic to expect them. When compared with results from inside r25, 
we find outer disk star formation to be distinct in several ways: it is extremely inefficient (depletion 
times of many Hubble times which are also long compared to either the free fall or orbital timescale) 
with column densities and SFRs lower than found anywhere inside the optical disks. It appears that 
the H I column is one of, perhaps even the key environmental factor in setting the SFR in outer galaxy 
disks. 
Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies — stars: formation 



> 

oo 

m 

o 
o 



% 



1. INTRODUCTION 

When galaxies are observed with sufficient sensitivity, 
star formation is often seen to extend well beyond the 
optical disks, reaching far into the extended H I disks. 
After early indication s of extended UV emission in 
nearby galaxies (e.g., iDonas et al.l |1981[ ). the discov- 
ery of extended UV (XUV) disks in a large number of 
nearby galaxies wa s one of the m ajor achieve ments of the 
GALE X mission (iThilker et al.l 2005; Gil de Paz et all 
2005"; 'Thilke r et al.l 120071: iGi l do Paz et al.' ■2007b; 



Boissier et al. l2007t IZaritskv fc Chr istlein 2007; 
Thilker et all 120091: faunter et al.l l2010l ). Extended 
star formation is seen not only in the UV, but also 
in the optical, e.g., in deep Ha and broad brand 
observations, which reveal populations of young stars 
in th e outer dis ks (Ferguson et al. 1998; Lclicvrc & Roy 
l2000t ICuillandre et al.l 120011 : Ide Blok fc Waited 120031: 
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Christlein fc Zaritskvl [20081: iHerbert-Fort et al.l I20T0I: 
Goddard et al.l l2010i: iWerk et al.l I2010D . This agrees 
well with the observations that indirect tracer s of past 
or li k ely future star f ormat i on, like du s t (iZaritskvl 
[T99I IPopescu fc T^ '2003'; 'Dong et al.l I2008h CO 
emission (iBraine fc HcrDin 2004; Brai ne et al.l 120071: 



Cardan et al. 
2007bt 



I2007f). and metals 

iBresolin et all |2009| ) 



are 



(|Gil de Paz et al.l 
detected in the 
extended H I envelopes of galaxies. 

Star formation at large radii usually does not account 
for a large fraction of a galaxy's total star formation rate 
(SFR), but studying this process offers a way to illumi- 
nate the ph ysics behind the star formation process (e.g., 
iBush et al.ll20Q8ll2010[ ): low metallicites and dust abun- 
dances, relatively high shear, low total gas column den- 
sities spread over significant scale heights, a preponder- 
ance of H I compared to H2 , and a comparatively weak 
stellar potential well make the ISM in outer galaxy disks 
a distinctly different environment compared to the typi- 
cal star-forming ISM in the inner part of a galaxy. This 
contrast leads to a much lower rate of star formation per 
unit gas mass at large galactocentric radii and to the 
conclusion that the physics behind the conversion of gas 
into staisjnust be affected by these environmental factors 
(e.g.,: Leroy et al. 2008). 

A robust, quantitative picture of how the environ- 
ment in outer disks affects star formation is important 
if we want to understand the origins of galaxy struc- 
ture. Star formation at large galactocentric radii will 
affect how chemical enrichment varies across a galaxy 
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(|Gil de Paz et al.ll2007b[ ) and plays a critical role in de- 
termining the location and form of the break in the ex- 
ponential stehar disk (iPohlen fc Truiillol [20061 ) . Also, 
many galaxies sustain a large reservoir of (low column 
density) gas in their outer disks over evolutionary time 
scales. Measuring the gas consumption time scale in this 
regime for many galaxies and comparing it to the gas 
consumption time found for the inner parts of galax- 
ies may provide valuable clues regarding the role of 
outer disk gas for fueling star formation over cosmologi- 
cal times (e.g., Shlosman. Fran k & Bcgclman 1989; Blitz 
[19961: iBauermeister et ani2009( ). 

In this paper, we study the relationship between 
atomic gas (H I) and star formation at large galacto- 
centric radii. We use state-of-th e-art HI ('The H I 
Nearby Galaxy Survey', THINGS, IWalter et al.[ MM) 
and UV data ('G ALEX Nearby Galaxy Survey', NGS, 
IGil de Paz et al.ll2007a) . which provide the field-of-view 
and the sensitivity needed to probe into the outer disks of 
galaxies while still offering the resolution to examine the 
interplay between gas (H I) and star formation locally, 
i.e., on scales of a few times 100 pc. 

We use these two datasets to study the relationship be- 
tween H I and star formation in 22 outer galaxy disks (de- 
fined as r = l-2xr25). We assess whether the observed 
decline of the SFR with galactocentric radius is predom- 
inantly due to a decreasing gas supply and we look for 
signs of star formation thresholds, as suggested theoret- 
ically and by observations of sharp truncations in radial 
distributions of H II regions (e.g., Martin fc Kcnnicuttl 
1200 It) . We examine radial and local variations in the 
star formation efficiency (SFE), i.e., the SFR normal- 
ized to the H I column (and thus the inverse of the gas 
consumption time) and variations in the spatial corre- 
lation of star formation and H I with galactocentric ra- 
dius. We compare our results for the outer disks of spirals 
to dwarf galaxies (which share many of the same envi- 
ronmental factors) and link them to observations from 
within the optical disks of an overlapping set of galaxies 
(jBigiel et al.. ,200a: ,Lcroy et al. 20^). 

2. DATA 

We study 22 galaxies: 17 spiral and 5 dwarf galaxies. 
This sample is constructed from the overlap of THINGS 
|Wahcr et al. 200 8) and targets of the GALEX NGS 
([Gil de Paz et al.l I2007al ) that were observed to similar 



depth (integration time of > 1.5 ks, corresponding to at 
least one orbit, which is the standard integration time 
for GALEX NGS targets). Table[T]lists our sample along 
with adopted distance, inclination, p osition angle, opti- 
cal radius r25 and morphology f from [Walter et aL[|2008l 
except that we adopt i = 20° in NGC 5194). We correct 
all maps for inclination using the angles given in Table 

m 

To allow a rigorous comparison, we degrade all H I and 
UV maps to a common resolution of 15" (set by the H I 
map with the lowest resolution) by convolving with a 
circular Gaussian beam. We have carried out a paral- 
lel analysis at a matched physical resolution of 1 kpc (set 
by the largest physical resolution in our sample) and find 
that for the spirals our results for the two cases differ only 
marginally. For the five dwarfs in our sample, working at 
1 kpc resolution significantly reduces the number of inde- 
pendent measurements (per galaxy and total) , constrain- 



ing our ability to robustly compare the two approaches 
due to the limited statistics. In the following we thus 
only present the matched angular resolution case. The 
average physical resolution in our sample is ~ 600 pc, 
varying between ^ 230 pc and ^ Ikpc across the sam- 
ple. 

We focus on the outer disks of galaxies, which we de- 
fined to be between 1 and 2xr25, where r25 is the isopho- 
tal radius corresponding to 25 B-band magnitudes per 
square arcsecond. This regime is illustrated in Figure 
[T] for one galaxy, NGC 3621. The left panel shows the 
THINGS H I and the right panel the GALEX far UV 
map. The overplotted annuli indicate 1 and 2xr25 in the 
plane of the galaxy, respectively. 

2.1. A Note on Terminology 

Throughout this paper we compare 21-cm intensity to 
far UV intensity. We will refer to 21-cm intensity and the 
surface densities of both H I and 'gas' interchangeably 
because we consider it likely that opacity effects in the 
21-cm line and the contribution of molecular gas are both 
small in the regime we study. Similarly, we measure far 
UV intensities, Ifuv, from the GALEX maps. We con- 
vert these intensities into an approximate SFR surface 
density (Ssfr) using the assumptions discussed below 
(§ 12. 4|) and will use the terms 'UV intensity' and 'star 
formation rate surface density' to mean the same thing 
(/puv)- In both cases, the reader interested in linking 
our plots directly to observables has only to make a lin- 
ear transformation of any axis using the equations given 
in this paper. 

When we refer to the conversion of H I into stars we 
implicitly assume that the gas becomes molecular first, 
i.e., H I forms into H2 which then forms stars. Because 
H2 in outer disks is not readily observable, we are forced 
to consider a 'zoomed-out' version of this process, the 
conversion of H I into stars with H2 as an unconstrained 
intermediate phase. H2 likely constitutes only a small 
fraction of the gas mass in outer disks, so the question 
of what drives the ISM to form stars in outer disks may 
still be robustly addressed using only H I and FUV. 

There is no set of wide-field CO maps t hat extends 
to_2r25 (the widest- field maps reach ~ r25 , iLerov et"all 
|2009( ). However, we can readily see that from the (aver- 
aged) star formation rate surface densities that we infer 
for outer disks (Ssfr « IQ-^-IQ-^Mq yr"! kpc"^ see 
§[3]) that CO is likely to be very faint and H2 is only 
a relatively minor part of the ISM. If these Esfr were 
found in the disk of a spiral galaxy, the corresponding 
H2 surface densiti es would be I]h2 '~ 0.01-1 Mq pc~^ 
(jBigiel et al.ll2008[ ). roughly corresponding to a single gi- 
ant molecular cloud (Mh2 ~ 10'*-10^ M©) per resolution 
element. This is low enough to assume that H I com- 
fortably dominates the ISM across our data on (approxi- 
mately) kpc scales and implies CO intensities well below 
the detection limits of most existing maps. Because of 
this negligible contribution of H2 to the gas budget and 
the lack of CO observations, we are confident using H I 
emission to trace the bulk distribution of the mass in the 
ISM in outer disks on kpc scales. 

Readers interested in comparing this paper to other 
results should also note that we neglect any contribution 
from helium or heavier elements when quoting gas sur- 
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TABLE 1 
Sample Properties-'- 



Galaxy 



D 

[Mpc] 



i PA 

[dcg] [deg] 



'■25 
farcminl 



^25 
[kpc] 



Hubble 
type 



Dwarfs 



DD0154 


4.3 


66 


230 


1.0 


1.2 


Irr 


Hoi 


3.8 


12 


50 


1.7 


1.8 


Irr 


HoH 


3.4 


41 


177 


3.3 


3.3 


Irr 


IC2574 


4.0 


53 


56 


6.4 


7.5 


SABm 


NGC 2366 


3.4 


64 


40 


2.2 


2.2 


Irr 


Spirals 


NGC 0628 


7.3 


( 


20 


4.9 


10.4 


Sc 


NGC 0925 


9.2 


66 


287 


5.4 


14.3 


Sed 


NGC 2403 


3.2 


63 


124 


7.9 


7.4 


SABc 


NGC 2841 


14.1 


74 


153 


3.5 


14.2 


Sb 


NGC 2903 


8.9 


65 


204 


5.9 


15.2 


SABb 


NGC 3198 


13.8 


72 


215 


3.2 


13.0 


Sc 


NGC 3351 


10.1 


41 


192 


3.6 


10.6 


Sb 


NGC 3521 


10.7 


73 


340 


4.2 


12.9 


SABb 


NGC 3621 


6.6 


65 


345 


4.9 


9.4 


SBcd 


NGC 3627 


9.3 


62 


173 


5.1 


13.8 


SABb 


NGC 4736 


4.7 


41 


296 


3.9 


5.3 


Sab 


NGC 5055 


10.1 


59 


102 


5.9 


17.3 


Sbc 


NGC 5194 


8.0 


20 


172 


3.9 


9.0 


Sbc 


NGC 5236 


4.5 


24 


225 


7.7 


10.1 


Sc 


NGC 5457 


7.4 


18 


39 


12.0 


25.8 


SABc 


NGC 7331 


14.7 


76 


168 


4.6 


19.5 


Sbc 


NGC 7793 


3.9 


50 


290 


5.2 


5.9 


Sed 



^ See [Walter et all 1200811 for further information on individual 
galaxies and for references to the values quoted in this table. 

face densities, but that we do take them into account 
when quoting star formation efficiencies or gas depletion 
times. 

2.2. THINGS HI 

To estimate the surface density of neutral, atomic hy- 
drogen, Ehi, we use VLA maps of the 21-cm line obtained 
as part of 'The H I Nearby Galaxy Survey' (THINGS, 
I Walter et aL|[2008l ). THINGS consists of high-resolution, 
high-sensitivity H I data for 34 nearby galaxies obtained 
with the NRACO VLA. We use 'natural' weighted maps, 
which offer the best possible signal-to-noise ratios and 
have an average synthesized beam size of ^ 11". These 
maps are sensitive (~ 3ct) to column densities of Shi ^ 
0.4 M0 pc~^ (relatively uniformly across our sample) 
at our working resolution of 15". Because the THINGS 
data include observations in the VLA's most compact 
(D) configuration, missing fl ux is not expected t o be a 
large concern. For details see lWalter et al.l ()2008D . 

The FWHM of the VLA primary beam (field-of-view) 
is 32'. In the cases of NGC 5236 (M83) and NGC 5457 
(MlOl), this limits the radii that we can consider to less 
than 2xr25; in these galaxies the analysis is carried out 
within 1.7xr25 and 1.5xr25, respectively. 

2.3. GALEX UV 

To trace recent star formation in the outer disks, we 
use far UV data from the 'GALEX Nearby Galaxy Sur- 
vey' (NGS, [GiIiilPaiiFinilOOZi). GALEX provides 
simultaneous imaging in a far UV (FUV, 1350- 1750 A) 
and a near UV (NUV, 1750 -2800 A) broadband filter 

® The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



with angular resolutions (FWHM) of 4.0" and 5.6" re- 
spectively, and a fiel d-of-view of 1.25° (for details see 
iMorrissev et al.ll2005D . We use only the FUV band to 
trace recent star formation, because it is less sensitive 
to any old stellar population and suffers less contamina- 
tion from foreground stars and background galaxies. We 
use the NUV band to identify foreground stars and to 
estimate the effects of extinction. 

We process the NGS maps slightly. First, we identify 
foreground stars as regions with an NUV/FUV inten- 
sity ratio > 10 and a signal-to-noise ratio > 5 in the 
NUV maps. We blank these regions after checking the 
results by eye. In a few cases we adjust the color cut-off 
to a value higher than 10 to remove particularly bright 
foreground stars. We also estimate and remove a small 
background from the FUV maps, which we measure away 
from the galaxy after discarding emission with intensities 
more than 3cr above the median value of the image. Fi- 
nally, we correct the maps for the effects of Galactic ex- 
tinction, which we estimate from the extinction maps of 
iSchlegelet all(fl 998) assuming Apuv = 8.24 x E{B-V) 
(|Wvder et &\Mmi\) . In NGC 5194 (M51), we also man- 
ually blank the interacting northern companion, M51b. 

A typical value for the RMS scatter of the noise in 
our FUV maps at our working resolution is ~ 2 x 
10^^ mJy arcsec^^ (l"')- The exact value varies signif- 
icantly from map to map and depends rather sensitively 
on the method used to derive it. This estimate refiects 
the median-based scatter (which avoids foreground stars 
and background galaxies) about the zero level in maps 
that have already been convolved to 15" resolution. Vari- 
ations in exposure time and map quality lead to a fac- 
tor of 2-3 scatter in the sensitivity of individual galaxies 
about the value quoted above. 

2.4. Conversion to Physical Units 

We convert 21-cm intensity into atomic gas surface 
density via: 



Em [Mq pc-2] = 0.015 /2icm [K km s^^] 



(1) 



which yields a hydrogen mass surface density and does 
not include heavy elements. To convert FUV intensi- 
ties in to SFR surface densities, we adopt the conversion 
of Sal im et al.l ()2007l thei r equation 10). For compati- 
bility with previous work (jLerov et al.ll2008l : IBigiel et alj 
|2008[) . we divide their coefficient by an extra factor of 
1.59, making the formula a pprop riate for a Kroupa-type 
IMF (compare iLerov et al.ll2008l Appendix D). Then 



SsFR [Mq yr-1 kpc"'] = 0.68 



X 10"^* X (2) 

^Fuv [ergs s~^ Hz^^ kpc"^] . 

iSaUm et all ()2007[ ) derived this calibration for --50,000 
galaxies by fitting population synthesis models to 
GALEX UV and SDSS optical multiband photometry. 
We assume here that the same calibration is applicable 
to the FUV intensity averaged over large parts of outer 
galaxy disks. A dire ct check on t his as sumption is not 
straightforward, but ILerov et"al] ()2008D compared vari- 
ous SFR tracers in the inner parts of many of our target 
galaxies and found a scatter of ~ 50%, which we will 
take as our uncertainty in Equation [2j We derive Esfr 
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Fig. 1.— THINGS Hi (left, linear stretch) and GALEX far UV (right, log stretch) maps for NGC 3621. The ellipses indicate 1 and 2 
optical radii (r25) in the plane of the galaxy. In this paper, we compare H I and far UV emission in this regime for 22 nearby galaxies. 



from the GALEX FUV maps alone, neglecting the ef- 
fects of internal extinction. This is both a practical and 
physical decision. Standard methods to correct for the 
effects of internal extinction — e.g., combination with IR 
intensity, use of NUV/FUV color — are largely impossi- 
ble because of the low signal-to-noise and easy confusion 
with background sources in this regime. At the same 
time, we expect the effects of extinction to be small, so 
that introducing a (noisy) correction factor could easily 
do more harm than good. 

We can estimate a likely upper limit to the bias intro- 
duced by neglecting extinction from our H I maps. In 
the Galaxy, N(R)/E (B -V) = 5.8 x lO^i cm^^ mag^i 
(jBohlin et al.l Il978| ). This should represent an upper 
limit to the extinction for a given H I column in outer 
disks, which have lower metallicities and dust abun- 
dances than the Milky Way. If we assume that FUV 
originates from the midplane (and is thus only af- 
fected by half the dust alon g a line-of-sight) an d that 
Apuv/E{B - U) = 8.24 (I Wvder et alT 120071) . then 
Apuv = N{HI) X 7.1 X 10-22 and the FUV intensity 
corrected for internal extinction is 



h 



FUV, con- 



h 



FUV 



^ -[^gJV(H/)x2.8xlO' 



(3) 



For an H I column of 3.0 Mq pc ^ (i.e. 



10 



3.7 X 
2°cm^2) ^ a typical value between 1 and 2 r25. Equa- 
tion [3] yields a correction factor of ^ 1.3. Thus we ex- 
pect an upper limit to extinction effects to be ^ 30%, and 
lower at larger radii, where H I surface densities drop sig- 
nificantly below 3.0Mqpc-2 and for any region with a 
dust-to-gas ratio lower than the Galactic average (almost 
certainly all of our disks) . 

While we expect that extinction is not a large concern 
in the outer disks, several plots will display data from 
within the optical disks (see, e.g., the radial profiles in 
i) l3.ip . In this regime, Esfr, as defined here, will repre- 
sent a (significant) underestimate of the true SFR surface 
density (due to dust extinction). As emphasized above, 
the conversion between Ssfr and the (observed) FUV 



intensity /puv is linear and it is thus straightforward to 
link any of the results or plots in this paper directly to 
our observables. 

3. RESULTS 

We employ several approaches to assess the gas-SFR 
relationship in outer disks. First, we compare radial pro- 
files of H I and FUV between 1 — 2 x r25 to test how the 
overall decline in gas and star formation relate to one 
another (ii l3.ip . Then, we measure the strength of the 
local correlation between EgpR and Shi as a function of 
galactocentric radius (§ 13. 2p . We subsequently investi- 
gate the gas-SFR relationship directly, testing how Esfr 
and the star formation efficiency depend on Shi in vari- 
ous regimes (i) l3.3|) and carry out a pixel-by-pixel analysis 
of the outer disk (ii l3.4l) . We assess the role of large scale 
gravitational instability as a potential driver for outer 
disk star formation (§ 13. 5|) and compare our outer disk 
data to resu lts for the optical disks of galaxies (§ 13. 6|) 
presented bv iBigiel et aD ()2008[ ). 

3.1. Radial Profiles and Exponential Scalelengths 

Both H I and FUV emission decline as a function of 
radius across the outer disks. Because gas is the fuel for 
star formation, one would expect the SFR to drop in a 
similar fashion compared to the H I (if all other factors 
were equal). To test this expectation, we derive radial 
profiles of each quantity and characterize the decline at 
large radii using an exponential scale length (i.e., the 
length over which the emission declines by a factor e). 

Figure [2] shows the radial profiles of H I (black) and 
FUV (gray) together, along with the exponential fits, for 
our 22 sample galaxies. The profiles are azimuthal av- 
erages over beam- wide (15"), concentric elliptical annuli 
constructed using inclination and position angles from 
Table [T] The error bar on each profile point indicates 
the RMS uncertainty in the mean in that ring. For easy 
comparison of profiles in different galaxies, we use com- 
mon Ehi and Ssfr axes. The vertical dotted line indi- 
cates the innermost data point included in the fit (which 
is the profile point closest to r25). 



Extremely Inefficient Star Formation in the Outer Disks of Nearby Galaxies 



Galactocentric Radius [kpc] 
0.0 0.5 1.0 1.5 2.0 



Galactocentric Radius [l<pc] 
12 3 4 5 5 



10^ 

10' 

10° 

10-' 



10" 



DD0154 



HI 
FUV 




- 10" 



- 10" 



-= 10" 



u 

Q_ 


O 
Q_ 







>^ 


^ 





I 



10' 



10' 



10° 



10" 



10" 



10" 



ilO" 



Holl 



0.0 0.5 1.0 1.5 2.0 

Galactocentric Radius [rjj] 

Goloctocentric Radius [kpc] 
1 2 




0.0 0.5 1.0 1.5 

Galactocentric Radius [rj^] 

Galactocentric Radius [kpc] 

" 10 12 14 




0.0 0.5 1.0 1.5 2.0 

Galactocentric Radius [rjj] 

Galactocentric Radius [kpc] 
12 3 4 



0.0 0.5 1.0 1.5 2,0 

Galactocentric Radius [''25] 

Galactocentric Radius [kpc] 
2 4 6 8 10 12 14 




0.0 0.5 1.0 1.5 2.0 

Galactocentric Radius [rjj] 



0.0 0.5 1.0 1.5 2.0 

Galactocentric Radius [''25] 



Fig. 2. — Radial profiles of SgFR (FUV, gray) and Spji (black) for the 22 galaxies in our sample. Error bars on each point show the 
(Icr) uncertainty in the mean in that ring and horizontal lines (black: Hi, gray: FUV) show conservative sensitivity estimates based on the 
linc-of-sight sensitivities in our maps and assuming 20 independent measurements per annulus (note that the FUV sensitivity is sometimes 
below the lower plot limit). A vertical line indicates the innermost data point included in the fit (approximately at r25). Solid lines show 
exponential fits to the decline between 1 — 2 X r25, only considering points above our sensitivity cuts. 
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Fig. [2] — continued. 
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Horizontal lines indicate the sensitivity limits: 
O.O9M0pc~^ for Shi (black), a value seldom reached 
within 2 x r25, and the sensitivity (2 a) of each FUV map, 
converted to units of Ssfr (gray). To estimate these sen- 
sitivity limits we use the H I sensitivity from i) l2.2l and the 
FUV sensitivity from each map and assume (a conserva- 
tive number of) 20 independent measurements contribut- 
ing to each annulus to reflect the gain in sensitivity by 
azimuthal averaging. 

Solid lines show our exponential fits to Ssfr (gray) 



and Shi (black) . We carry out the fits between ~ 1 (the 
vertical dotted line) and 2 x r25 and consider only pro- 
file points above our sensitivity cuts. We reiterate that 
for NGC 5457 and NGC 5236 the field-of-view of our 
H I maps is too small to probe out to 2 x r25 . We thus 
blank their radial profiles outside their respective limit- 
ing radius (compare ii l2.2p . We note that NGC 3627 has 
too few data points above the H I sensitivity cut to fit 
the profile beyond ^25. We also note that in particular 
NGC 5236 (M83) is known to have extended FUV fea- 
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tures in the outer disk (e.g.. iThilker et ai1l2005D . These 
features have a very low fiUing fraction for the outer disk 
annuh and are thus not easily visible in our radial pro- 
files. This galaxy is subject of a separate paper where 
significantly deeper and more extended data are used to 
assess the relationship be tween H I and FU V emission 
out to many optical radii (jBigiel et al.ll2010() . 

To fit the profiles we use an ordinary least-squares 
(OLS) approach and estimate the uncertainties in the 
resulting fit in several ways: by adding noise and re- 
fitting, bootstrapping (resampling with repeats), vary- 
ing the radii used to define the fit, and considering only 
a subset of azimuthal angles (rather than the whole ring) . 

In fact, azimuthal variations within a galaxy appear to 
dominate the uncertainty in the fit, which might be ex- 
pected from the common appearance of streamers, tidal 
features, and outer arms in both the H I and SFR maps. 
To separate out this effect, we derive four independent 
scale lengths for each galaxy, one for each of four sec- 
tors with a 90 degree opening angle. The remaining fac- 
tors (bootstrapping, radius definition, noise) contribute 
to an uncertainty in the scale lengths for each sector. 
The typical scatter among the four sector scale lengths 
is ~ 0.1 X r25, with significantly less scatter in the FUV 
profiles than the H I (indicating a higher degree of az- 
imuthal symmetry in the FUV emission than in the H I) . 

Figure [3] shows H I scale length as a function of FUV 



1 0"^ . , , , I , , , , i , , , , I , ' , : ' ■ 1 0" 
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scale length. Each galaxy contributes four points, one 
from each 90°-wide sector, with points for a galaxy shar- 
ing the same color and symbol. Black points represent 
dwarf galaxies, otherwise color is arbitrary. Four galax- 
ies (DDO 154, NGC 2841, NGC 3198, NbC 3621) and 
9 individual sectors are omitted from this plot because 
their H I profiles are too fiat to be fit robustly over the 
range we consider. Their profiles do appear to decline 
roughly exponentially, but the derived scale length is 
> 1 X r25 , yielding too small a dynamic range for a robust 
fit . NGC 3627 is also omitted as there are too few data 
points above the sensitivity limit to fit the H I profile (see 
above) . 

Figure [3] demonstrates that in almost every outer disk 
quadrant we study, the surface density of star formation 
traced by FUV emission drops more quickly than the 
surface density of gas (for comparison, the dashed line 
shows equality). An exception is one sector in NGC 925 
(filled purple diamond) which shows a slightly larger 
FUV scale length (though with large associated uncer- 
tainties). Even though generally there is a steady de- 
cline in both H I and FUV with radius (Figure ^ and 
a clear correlation between the scale lengths of the two 
declines, these declines are not identical. Instead star 
formation drops with roughly half the scale length of the 
H I. This is illustrated by the solid line, which shows 
the mean ratio of FUV scale length and H I scale length 
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Fig. 3. — Scale lengths for the exponential decline of FUV (y-axis) and Hi (x-axis) in outer galaxy disks. Each galaxy contributes four 
independent measurements from different azimuthal sectors; data points that belong to the same galaxy share the same color and symbol. 
Black data points represent dwarf galaxies. The dashed line indicates a slope of unity. Virtually all data lie underneath this line, indicating 
that the FUV profiles generally decline more rapidly than the Hi profiles. The solid line, which shows the mean ratio of the two scale 
lengths, has a slope of ~ 0.5, i.e., the outer disk Hi scale length is twice as large as the FUV scale length on average. 



of 0.5 (the la scatter about the mean is ^ 0.2). There- 
fore, and despite the remarkable extent of star formation 
revealed by GALEX, it is worth bearing in mind that 
this widespread UV emission is still relatively centrally 
concentrated compared to the H I. 

3.2. Local Correlation 

Comparing radial profiles tells us about the bulk be- 
havior of H I and FUV in outer disks, but does not reveal 
information about any local relationship. To assess the 
detailed relationship between FUV and H I, we now turn 
our attention to how the two quantities correlate on a 
line-of-sight by line-of-sight basis at a given galactocen- 
tric radius. Our approach is to divide each galaxy into 
10 equally wide radial bins between the galaxy center 
and 2 x r25 and then compute the Spearman rank cor- 
relation coefficient between FU V and H I in eac h bin. 
The rank correlation coefficient (jPress et al.lll992n . r, is 
a non-parameteric measure of the strength of any one-to- 
one relationship between two quantities. Possible values 
range from r = —1 to 1, with r = 1 indicating a perfect 
correlation — i.e., that the brightest data point in H I is 
associated with the brightest data point in FUV, the sec- 
ond brightest in H I with the second brightest in FUV, 
etc. — but giving no information on the functional form. 
On the other hand r = indicates a lack of correlation 
(expected if the two distributions are independent), while 
a perfect anti-correlation will yield r = — 1 . 

Figure |4] shows r as a function of normalized galacto- 
centric radius for the spiral (left) and the dwarf (right) 
galaxies separately. Each gray data point represents r 
measured for one galaxy in one 0.2 x r25-wide radial bin 
(for better visibility we add a small amount of noise to 
the x-position to distinguish individual measurements in 



the same bin from one another). The vertical error bars 
show the scatter in r derived from randomly pairing H I 
and FUV data for that bin; this process destroys any 
correlation, yielding r w and the scatter observed re- 
peating the process 100 times gives us an estimate of the 
uncertainty (this is most rigorously thought of as how 
confident we are that r j^ 0). The black line condenses 
the gray points into a single trend by giving the median 
correlation coefficient for all galaxies in each radial bin 
with error bars indicating the la RMS scatter among 
galaxies. 

For spiral galaxies, r rises from ~ 0.1 near galaxy cen- 
ters to ~ 0.7 at about r25 and then drops again to ~ 0.3 
near 2 x r25- To assess whether the low values of r at 
very small and very large radii still represent statistically 
significant correlations, we extend the procedure that we 
used to estimate the errors for the individual measure- 
ments. We randomly pair H I and FUV data (now for 
all spirals and dwarfs, respectively), so that r = by 
construction, and bin these data and then derive a from 
the scatter in r about the known true value of 0. This 
exercise leads us to estimate la < 0.04 for the median 
profiles in both the spiral and dwarf sample, implying 
that the median r is significantly different from in the 
outer parts of spirals but is consistent with within the 
uncertainties in the inner parts of spirals and the outer 
parts of dwarf galaxies. 

Inside '^ 0.2 x r25, the distributions of H I and FUV are 
compatible with the hypothesis that r — 0, i.e., that the 
two quantities are unrelated. This regime is not the focus 
of this study, but the result is easy to understand given 
that we neglect internal extinction, which is significant 
in this regime, and molecular gas, which dominates the 
gas budget in the ISM at these radii. 
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Fig. 4. — Spearman rank correlation coefficient, r, relating Hi and FUV as a function of galactocentric radius in the spiral (left) and 
dwarf (right) galaxies. For each galaxy, we measure r in 10 radial bins spaced evenly between the galaxy center and 2xr25. Each of these 
measurements, r for one radial bin in one galaxy, appears as a gray point. Vertical error bars show Monte Carlo uncertainty estimates. 
The black curve in each panel indicates the median r over all galaxies in each radial bin. Black error bars show the Icr RMS scatter among 
galaxies in that bin. FUV and Hi emission are significantly correlated, with r > 0.5, around the optical radius in the spirals and throughout 
the optical disks of the dwarf galaxies. In the inner parts of the spirals extinction and the (unaccounted for) presence of H2 weaken the 
correlation. Outside the optical radius, r decreases with increasing radius in both dwarfs and spirals. 

The effects of molecular gas and internal extinction 
should be negligible in outer disks. The decrease in r 
from its peak near the optical radius towards larger galac- 
tocentric radii appears to represent an intrinsic weaken- 
ing of the correlation between star formation and total 
gas in this regime. 

Maps of both FUV and H I in outer galaxies are often 
dominated by tidal streamers, arms or other relatively 
narrow but extended features. One way that the result 
in Figure |4] might arise even in the face of a one-to-one 
correlation between star formation and gas is if such fea- 
tures exist but are systematically offset between the H I 
and FUV maps. To check this possibility, we repeat the 
analysis in this section with a set of maps convolved to 
30" resolution. If a small systematic offset were driving 
the behavior at large radii in Figure IH we would expect 
r to increase at lower resolution. We do not find such 
an effect, suggesting that systematic offsets in otherwise 
perfectly matched maps do not drive the observed decline 
in r with radius. 

Dwarf galaxies (right panel, Figure HJ do not exhibit 
the depression in r at small radii seen in spirals, which 
can be explained because extinction and the contribution 
of H2 to the total gas budget are generally small even in 
the inner disks of dwarf galaxies. It is not clear if the 
slightly lower r at large radii found in dwarfs compared 
to spirals is real or a statistical artifact: our sample of 
dwarfs is small (5) compared to our sample of spirals (17) 
and the galaxies themselves are smaller, leading to fewer 
independent line-of-sight measurements per galaxy. The 
outer disks of dwarf galaxies do show, however, a decline 
in r with increasing radius similar to that seen in the 
spirals. 

3.3. The Star Formation Efficiency in Outer Disks 

We have seen that H I and FUV emission exhibit dif- 
ferent radial behaviors (Section I3.ip and that the local 



correlation between the two declines from good (r w 0.7) 
near r25 to poor (r w 0.3) near 2 x r25. The first result 
means that the rate of star formation per unit gas changes 
systematically across outer galaxy disks while the second 
implies that quantities other than the local gas reservoir 
may be important to set this quantity. Here we directly 
investigate the variation of the rate of star formation 
per unit gas with H I column and radius in outer galaxy 
disks. Following convention for extragalactic studies, we 
refer to the star formation rate per unit gas (here FUV- 
per-H I) as the star formation efficiency (SFE). This is 
the inverse of the gas depletion time {t^^ ), the time 
required for present day star formation to consume the 
available gas reservoir; it is sometimes quoted as a true 
(dimensionless) efficiency by normalizing to 10® yr (i.e., 
quoting the fraction of gas consumed every 10® yr). The 
three formulations are all equivalent within a constant. 
AU SFEs (depletion times) quoted in this paper include 
a factor of 1.36 to account for heavy elements. 

In the remainder of this section, we analyze a data set 
of matched H I and FUV intensities, each measured over 
an individual line of sight at 15" resolution and together 
covering the area between 1 and 2 x r25 for our whole 
sample. We extract th ese data from the maps following 
the approach used by iBigiel et"aI1 (|2008l ) to study star 
formation inside the optical disk. Briefly, we measure 
the intensity from the maps (convolved to 15" resolution) 
at non-overlapping sampling points (i.e., separated by 
approximately a beam width) and spaced evenly to cover 
the area between r25 and 2 x r25 in each galaxy. Each 
data/sampling point is assigned a weight equal to the 
inverse of the total number of data points for the galaxy 
from which it is drawn. We apply these weights when 
combining data in order to give equal weight to each 
galaxy; this avoids a few large ga laxies domin a ting a ll of 
our plots. Other details are as in IBigiel et al.l (J2008D . 
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3.3.1. Distribution o/Hl Columns 

Before we examine how the SFE depends on H I col- 
umn density, it will be useful to see the actual distri- 
bution of H I columns in our data. We show this in 
Figure [5] via normalized Shi histograms for the spiral 
(left panel) and the dwarf (right panel) samples. The 
dotted histograms show the relative distribution of H I 
columns across the entire outer disks (between 1— 2xr25), 
where we divided the H I surface density range from 
0.4 — 15 M0 pc~^ into 10 equally wide bins (there are 
very few data points at higher columns and the lower 
limit corresponds to the sensitivity of the H I maps). In 
constructing the histograms, we give equal weight to each 
galaxy rather than each data point. We test the robust- 
ness of the histograms by repeatedly adding noise to the 
original data and rebinning. Noise and false positives 
from regions below the sensitivity cut only scatter the 
results by a few percent. 

Unsurprisingly, the dotted histograms in Figure [5] 
show that low H I surface densities dominate outer 
galaxy disks. The 50*'' percentile surface density is 
~ 1.6 M0 pc~^ for spirals and ^ 2.3 M© pc^^ for dwarf 
galaxies. The spiral histogram appears approximately 
exponential, with an e-folding every ~ 1.6 M© pc~^. The 
dwarf histogram, which is based on fewer data, shows a 
shallower decline with a suggestion of a flattening at low 
columns. Broadly, higher column densities are more com- 
mon in the outer parts of dwarf galaxies than in spirals. 

We will also be interested in how the distribution of 
Ehi varies with radius. The other histograms in Figure[5] 
show this. We divide the data into three radial bins (0.5- 
1.0 r25, 1.0-1.5 r25, and 1.5-2.0 r25) and then separately 
plot the same kind of normalized histograms for each bin. 

The relationship between dwarfs and spirals does not 
change dramatically with radius. Spirals (left) show a 
steeper distribution of Shi than dwarfs (right) at all 
radii, so that dwarfs always have more high-column H I 
than spirals. Both spirals and dwarfs do show a signif- 
icant evolution in Ehi with radius: The fraction of low 
H I surface densities increases strongly with increasing 
galactocentric radius and at the lowest radii, Ehi flat- 
tens. The distribution even turns over for dwarfs, so 
that low H I columns are not the most common values 
(instead Ehi ~ 4 — 6 M© pc~^ is the most common value 
in dwarfs between 0.5 and 1.0 r25). 

3.3.2. SFE vs. H I Column 

Figure H] shows the SFE as a function of Ehi for data 
between 1 and 2 x r25. We bin the data by Ehi, assign 
equal weight to each galaxy and plot the median SFE 
and Ehi in each bin. Error bars indicate the Icr RMS 
scatter among the data. Very few data contribute to the 
highest column bins and as a result we cannot robustly 
estimate the scatter there. 

The measurements in Figure [B] appear to all be statis- 
tically significant. We check this by repeating the mea- 
surement after substituting a pure noise map for the FUV 
map and by directly propagating the statistical uncer- 
tainty in the two maps to derive the uncertainty in the 
median for each bin. A pure noise measurement would be 
well off the bottom of the plot and the formal uncertainty 
in each bin is extremely small (as might be expected from 
the large number of data). 



The figure shows that aside from a slight offset the 
outer disks of the dwarfs and the spirals exhibit a very 
similar SFE at a given H I surface density. In both cases, 
the SFE is a clear function of Ehi, increasing with rising 
H I column density. Power law fits to the points in Figure 
[6](solid lines) yield indices of 0.7±0.1 for both spirals and 
dwarfs. Although we saw above that the distribution of 
Ehi does differ slightly between the outer disks of spiral 
and dwarf galaxies. Figure [6] supports the idea that the 
outer disks of spirals and dwarf galaxies represent similar 
regimes regarding the regulation and timescales of star 
formation and that Ehi niay be key to this regulation in 
this regime. 

Stud ying a subset of our s piral sample, iLerov et"al] 
(|2008D and lBigiel erall (|2008D found a constant SFE for 
the H2 that dominates the gas reservoir in t he inner part 
of spiral galaxies. The value measured by iBigiel et al.l 
(l2008i') appears as a d ash-dotted horizontal line in Figure 
[51 ILerov et all ()2008[ ) further found the SFE of the total 
gas to be a well-defined function of galactocentric radius 
in both spiral and dwarf galaxies; we plot the SFE at 
'"25 predicted by their fit for each subsample as dashed 
lines. Comparing our current data to these lines we see 
that the SFE in outer disks is an order of magnitude 
or more below that of the star-forming gas in the inner 
parts of galaxies, so that it will take well in excess of a 
Hubble time to consum e the gas res ervoir in situ. The 
r25 predictions from the ILerov et al.l (2008) fits coincide 
with the high end of our measured SFE, suggesting that 
the steady decline they find may continue past r25. 

3.3.3. Dependence on Radius 

ILerov et al.l (|2008[) and IBigiel eFaLJ (|2008D found that 
galactocentric radius was actually a much better predic- 
tor of SFE than Egas in the inner parts of galaxies, likely 
because it tracks key environmental quantities like the 
stellar potential well and metallicity. To assess whether 
a similar radial gradient for the SFE holds in outer disks, 
we separate our data into three radial bins — 0.5 — 1 x r25, 
1 — 1.5 X r25 and 1.5 — 2 x r25 — and repeat the above 
analysis (? I3.3.2|) for each bin. 

Figure [7] shows the SFE as a function of Ehi for each 
radial regime with dwarfs and spirals plotted separately. 
All plot parameters are identical to Figure [HI Again we 
plot no error bars if only a few galaxies contribute to a 
bin. We also plot arrows instead of error bars where the 
scatter exceeds the lower plot boundary. Reading Figure 
[3 it is important to bear in mind that the SFE we plot 
between 0.5 - 1.0 r25 in spirals is biased (in opposite di- 
rections) by our failure to account for internal extinction 
or H2 (in these galaxies the ISM appears to be roughly 
equal parts H I and II2 at ^ 0.5 r25, Leroy et al. 2008); 
the points in Figure [7] are rigorously FUV-per-H I rather 
than SFE. 

In order to understand the upturn at low H I columns, 
we traced the data for the lowest H I-column bins back 
to the original images and found that they mostly arise 
from local H I depressions (i.e., 'H I holes') in regions 
with diffuse, locally smooth FUV emission. The FUV 
intensities involved are very low and inside a galaxy disk 
- the effect appears most pronounced for the inner radius 
bin in Figure [71 i.e. within r25 - interpreting these inten- 
sities can be complicated; one might expect some level 
of FUV emission from any galaxy due to an intermedi- 
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Fig. 5. — The 4 histograms in each panel (left: spirals; right: dwarfs) show the normalized Hi histograms between 1 and 2 X r25 (black 
dotted line) and for 3 other radial regimes: 0.5 — 1 X r25 (thick black line), 1 — 1.5 X r25 (black hashed) and 1.5 — 2 X r25 (filled gray). Every 
galaxy is assigned equal weight. In particular for the spirals, there is an approximately exponential decline in frequency with increasing 
Sjji, with dwarfs showing a much shallower decline and more high-column data. Each dwarf histogram is systematically shallower than 
the respective one for the spirals and with increasing radius, the fraction of sampling points with low Hi columns increases strongly at the 
expense of those with high Hi columns. 
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Fig. 6. — SFE (i.e., SFR/gas or FUV/Hl) binned by Shi for both samples (spirals: gray; dwarfs: black). Individual filled circles represent 
the median SFE and Sjji, error bars the Icr RMS scatter in each bin (error bars are omitted on the highest two bins because only little 
data contribute). Every galaxy is assigned equal weight. The solid lines indicate power law fits to the two distributions (both cases have 
slope r~ 0.7). Horizontal lines illustrate various comparison measurements: the SFE of molecular gas within the optical disks of spirals and 
of total (atomic plus molecular) gas at r25 in dwarfs and spirals. The SFE is lower in outer disks than within the optical disk and rises 
with increasing Hi column, showing a similar behavior in dwarfs and spirals. 
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ate age population, dust-scattered FUV light, or evolved 
stars. Therefore, while Figure [7] shows the real trend in 
our data we emphasize the systematic uncertainty in the 
lowest bins, particularly at low radii. 

At higher Shi (^ 5 M© pc~^) all radial regimes in both 
subsamples exhibit a positive correlation between SFE 
and Efji similar to what we saw in Figure [S] for the entire 
outer disks. The other general trend evident for spirals 
is a decrease in SFE with increasing radius at fixed Eri- 
For the dwarf sample the radial behavior is much less 
clear. A decline in FUV per H I with increasing radius 
is clearly present at low Shi, but whether this represents 
a real decline in SFE depends on one's interpretation of 
the low FUV intensities. 

3.4. Local Relations in Outer Disks 
3.4.1. Star Formation and H I 

We have so far binned and averaged our data to com- 
pare H I and FUV in a number of ways. It is also in- 
teresting to examine how individual data are distributed 
in EsFR-5]gas (FUV-H I) parameter space. Here we plot 
FUV as a function of H I directly for the data assembled 
in Section 13.31 This allows rea dy comparis o n to work 
on the inner parts of galaxies by iBigiel et"al] (|2008[ ) and 
numerous literature measure ments aimed at c onstraining 
the star formation law (e.g.. iKennicuttjflQOSf) . 

We have many thousands of data points, making a di- 
rect scatter plot impractical. Therefore, we convert the 
data into two-dimensional density distributions, again 
giving each galaxy equal weight. Figure [5] shows the re- 
sulting data density using contours to indicate the area 
containing the densest 25% (magenta), 50% (red), 75% 
(orange), and 90% (green) of the data. The data distri- 
butions in Figure |8] and the corresponding distributions 
inside r25 are available as supplemental online material. 

Many of the conclusions from Sections S I3.3.1I and 
ii l3.3.2l are again evident in Figure HI Depletion times 
are large (lines of constant H I depletion time appear as 
dotted diagonal lines in Figure [8]) and change systemat- 
ically but relatively weakly with changing Ehi. Dwarf 
galaxies exhibit somewhat higher Shi than spirals, lead- 
ing to a lack of low-column points in the right panel of 
Figure [S] At a given H I column density, the FUV one 
finds in spirals and dwarfs is quite similar. This last con- 
clusion can be clearly seen from the right panel of Figure 
[51 where the orange contour from the left panel appears 
as a thick black contour that closely matches the distri- 
bution observed in dwarfs. 

Sensitivity is a significant concern in this plot. The 
horizontal line shows a typical 3cr sensitivity for our FUV 
maps. A large fraction of our measurements lies below 
this line. This is problematic for a log-log plot, where 
negatives are not reflected. To robustly follow the gen- 
eral trend down to low Ssfr, we overplot median val- 
ues for EsFR in 5 equally spaced Shi bins as black cir- 
cles. All data, including negatives, contribute to the me- 
dian, making it much more sensitive than each individual 
point. Error bars on these points give our best estimate 
for the intrinsic (log) scatter in Ssfr in each H I bin. We 
derive this estimate by comparing the observed data in 
each bin to a series of mock data distributions. These 
are constructed to have the observed median and appro- 
priate Gaussian noise (measured from the FUV maps) 



with varying degrees of lognormal scatter (from 0.0 to 
2.0 dex). We compare each mock distribution to the ob- 
serve d data using a two -sided Kolmogorov-Smirnov (KS) 
test (iPress et al.lll992 1. The scatter that generates the 
distribution most similar to the observed data is our best 
estimate for the true intrinsic scatter. 

For H I columns > 3 M© pc~^, the scatter for both 
spirals and dwarfs is a few tenths of a dex (la), so that 
knowing the H I column (and that one was in an outer 
galaxy disk) could serve as an estimate for Ssfr with 
an accuracy of a factor of 2-3. While this continues to 
be true at low columns in dwarf galaxies, the scatter in- 
creases significantly below Shi ~ 3 Mq pc~^ for spirals 
(although the median Ssfr for spiral and dwarf galaxies 
is still roughly the same at these columns). If the ob- 
servation is real and not an artifact of poor statistics in 
the dwarf sample it may reflect a wider range of environ- 
ments in the outer disks of spiral galaxies. 

One feature that particularly stands out in Figure El 
but was not evident earlier, is an apparent maximum 
to the SFE in outer galaxy disks. Almost all data have 
H I depletion times longer than ~ 10^" yr (this line is 
a rough upper envelope to the green contour in both 
panels) . This implies that at their present SFR, even the 
most efficient outer disks in our sample will not consume 
their available gas in less than 10 Gyr (roughly a Hubble 
time) at their present SFR. 



3.4.2. SFE, Radius, and H I 

In Figure [21 we use the data to plot the SFE (t^^M as a 
function of radius. The corresponding plot showing SFE 
as a function of H I for individual lines of sight appears as 
Figure [ini In both plots contour levels and other details 
are identical to Figure [51 Black filled circles and error 
bars show the median SFE and best-estimate intrinsic 
scatter (using the same approach as above) in evenly 
spaced bins. 

Again several earlier conclusions are more clearly vis- 
ible in Figures M and [TUl the SFE (rocp) does appear 
to vary systematically with both radius and H I column, 
but this variation is quite weak compared to the intrin- 
sic scatter in the data (especially for the radial varia- 
tions in spirals). The difference between the dwarf and 
spiral samples seems to be largely the distribution of 
H I columns: at a given column or radius the SFE ap- 
pears largely similar in the two samples and the range of 
SFEs found for the two samples is similar (tdgp ~ 10^° 
- 10^^ yr), with 10^° yr a rough lower limit to Toop. 

These line-of-sight measurements also allow us to look 
for a star-formation threshold, a set of local conditions 
below which star formation is suppressed and above 
which it sets in. In a plot like in Figures [51 or [TUl a true 
threshold as a function of the ordinate should correspond 
to a rapid decline in the SFE as one crosses the thresh- 
old. Neither Figure [51 nor Figure [TUl show clear evidence 
for such a behavior; the SFE varies systematically but 
does not drop precipitously at any particular H I column 
density or radius. This does not necessarily rule out a 
local threshold, but if one exists then the timescales and 
spatial scales involved must be long enough and small 
enough to yield a smooth trend of FUV/H I as a function 
of H I column and radius (albeit one with large scatter). 
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Fig. 7. — SFE as a function of Sjji for spiral (left) and dwarf (right) galaxies. Methodology as Figure [6] except that the data have now 
been divided into 3 radial bins with results for each bin plotted separately (black filled circles show radii 0.5 — 1 X r25, dark gray circles 
from 1 — 1.5 X r25 and light gray circles from 1.5 — 2 X r25). We also plot arrows instead of error bars where the scatter exceeds the lower 
plot boundary. Generally, the SFE increases with Shi a-nd for a given Sjji, the SFE decreases with increasing galactocentric radius. 
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Fig. 8. — Pixel-by-pixel distribution of FUV (right axis; left axis after conversion to SgpR, Equation [2} as a function of Hi in the outer 
disks (1 — 2 X r25) of spiral (left) and dwarf (right) galaxies. Contours show the density of data after combining all galaxies in each sample 
with equal weight given to each galaxy. Magenta, red, orange, and green areas show the densest 25%, 50%, 75% and 90% of the data, 
respectively. Dotted lines indicate constant Hi depletion times of 10* to lO^'^ yr (taking into account heavy elements). A horizontal dashed 
line indicates the typical 3cr sensitivity of an individual FUV measurement. Black filled circles show our best estimate for the true relation 
between FUV and Hi after accounting for finite sensitivity: they represent the median FUV after binning the data by Shi and error bars 
are the lognormal scatter that yields the best match to the data after accounting for noise (see text). To allow easy comparison, we overplot 
the orange (75%) contour for the spirals as a thick black contour in the dwarf (right) plot. 
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Fig. 9. — FUV-to-Hl ratio (SEE or t^^ , y-axis) as a function of galactocentric radius in the outer disks of the spiral (left) and dwarf 
(right) galaxies. Contour levels, medians (black circles) and estimated intrinsic scatter (black error bars) are derived as for Figure [8l (also 
compare text). As for the SFR-Hl plots in Figure |8] both distributions are found to look almost identical, supporting the previous finding 
that the outer disk data from both galaxy samples lack clear distinguishing characteristics regarding their star forming properties. The 
SEE for both samples shows only mild variations across the outer disks. 
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Fig. 10. — The SEE {tj^^ ) as a function of Hi surface density Sjji for outer galaxy disks (1 — 2 X r25) in spiral (left) and dwarf (right) 

galaxies. Contour levels, the methodology used to derive the median SEEs (black circles) and associated scatter, and other plot parameters 
are identical to the previous figures. We do not find an obvious discontinuity in the the SEE as a function of Ejji (as one might expect for 
a star formation threshold). Instead the SEE changes relatively smoothly as a function of Ehi with significant scatter. 
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3.5. Toomre's Q in Outer Disks 

Toomre's Q parameter (jToomrdll964[) . which measures 
the stabihty of a thin axisymmetric disk (or ring, if the 
galactocentric radius is kept as a free parameter), has 
been closely linke d to the decline of star fo rmation in 
outer disks (e.g., iMartin fc KennicuttI |200TI ). Here we 
compare Q to the rate of star formation per unit gas in 
outer galaxies to see if a universal behavior emerges. We 
calculate Q via 



^ gas 



TrCrSnas 



(4) 



with values < 1 indicating instability (i.e., gas can col- 
lapse and form stars) and values > 1 indicating stabil- 
ity (i.e., no star formation). Here a gas denotes the ve- 
locity dispersion of the gas, G is the gravitational con- 
stant and K is the epicyclic frequency. In the outer 
disks that we study, a flat rotation curve is usually a 
good approximation fsee lde Blok et al.ll2008[ ) and we can 
calculate k from k = lAlvfiat/rgai, where Vd^t is the 
rotation ve locity at large radii. We adopt Uflat from 
iLerov et al . (2008) where available and derive it in a simi- 
lar manner for galaxies that they did not study. We take 
Cgas ~ llkm/s everywhere, which is a reasonable ap- 
proximation for the H I velocity dispers ion (to within a 
few k m/s) over outer galaxy disks (see iTamburro et al.l 
120091 ). Equation H] assumes the mass budget to be domi- 
nated by gas (ignoring the stellar contributio n considered 
by, e .g., Jog fc Solomon 1984: Rafikov 200lt lYang et al.l 
|2007() ; we consider this a reasonable approximation out- 
side 7-25. 

Figure [TT] shows the FUV-to-H I ratio as a function of 
Q. Low values of Q indicate gas that is more unstable 
to collapse and might thus be expected to correspond to 
regions of more active star form ation per unit gas (e.g., 
iLi et all [20051: lYang et al.ll2007D . Such a trend is not 
obvious from Figure [TTJ which shows that on the scales 
we study most gas appears quite stable with little or no 
local correlation between Q and the FUV-to-H I ratio. 

Despite the lack of a strong local trend, there is an 
overall correspondence: the depletion times that we de- 
rive are much larger than those found inside r25, while 
the typical values of Q in outer disks {Qmedian ~ 30) 
are, on the w hole, much larger t han those found inside 
r25- Inside r9.^. ILerov et al.l (|2008[ ) found Qmedian ~ 4 for 
the same methodology that we use here an d QmfAian. ^ 2 
considering a disk of ga s and stars (also see lBoissier et al.l 
l2003HYang et al.ll2007[ ). Figure HI] does not offer clear ev- 
idence for Q as a 'silver bullet' for star formation thresh- 
olds, but outer H I disks are clearly more Toomre stable 
than gas inside r25 (for a discussion of the interplay be- 
tween H I phases, the g as velocity dis persion, and Q see 
ISchavel[2004: dc Blok fc"WalteHl2006n . 

There are several subtleties to calculating Q: correc- 
tions are sometimes applied for disk thickness (stabiliz- 
ing) and the influence of stars (destabilizing) and the 
appropriate H I velocity dispersion is a matter of some 
dispute. We use the median velocity dispersion measured 
at r25, but using a lower dispersion more appropriate to 
a cold phase (e.g., agas = 6 km/s instead of llkm/s) 
would not change our basic conclusions, as Q would only 
decrease by a factor of '^ 2 — 4 (even before accounting 
for the fraction of gas in the cold phase, which must be 



relatively low given the observed dispersions). 

We do also note that there is a signiflcant contribution 
of dark matter to the local mass volume density in the 
outer disks. Assuming a flat rotation curve, the local 
dark matter volume density pdm at 1.5 r25 is of the same 
order as the local H I density, i.e., pdm ~ Phi- Deriving 
the effect of the dark matter on the disk stability is non- 
trivial, but even if it were concentrated into a thin disk 
along with the H I it would not be enough mass to induce 
formal instability. 

3.6. Comparison to Optical Disk Measurements 

We have looked at the relation between H I and FUV 
in outer disks in some detail. What remains now is to try 
to link these results to the inner parts of galaxies. We do 
so in the fo llowing by comparing our resu lts to those of 
ILerov et al.1 (|2008l ) and 'Bigiel ctal.' (2008), who studied 
an overlapping sample of spirals and dwarfs using many 
of the techniques used here. It is absolutely essential 
to note that as we do so, we abandon the clean map- 
ping between physical quantities and observables main- 
tained throughout the first part of this paper: Egas in 
ILerov et al.] (|2008f ) and Bigiel et al. (2008) includes the 
contribution of H? estimated from HERACLES CO maps 
(jLerov et al.ll2009l ) and Esfr represents a combination of 
FUV and IR (24/Lim) emission, with the IR usually the 
dominant term (see those two studies for details). With 
this caution firmly in place, we now compare the scale 
length of star formation, the local relation between Ssfr 
and Egas, and the SEE as a function of radius across the 
two regimes. 

Radial profiles of Ssfr appear to follow a steady de- 
cline with approximately constant scale length from in- 
ner to outer disks. We sh ow this in F igure [T2| where 
we plot scale lengths fit by ILerov et al.l (2008) to com- 
bined UV and IR profiles against our outer disk mea- 
surements. There is significant scatter in the plot, but 
with two notable exceptions the data appear to scatter 
around equality with ±0.1 r25 (typically ^ 50%) scatter. 
The two exceptions are the dwarf irregulars IC 2574 and 
Ho I, which both show a sharp downturn in their H I and 
FUV profiles starting just inside r25 (see Figure [2]). 

By contrast, just glancing at Figure [2] makes it clear 
that the H I does noticeably change its behavior between 
the inner and outer disk. Our H I profiles are usually 
fairly flat inside the optical disk (sometimes with cen- 
tral depressions), while almost all of the profiles show 
a radial decline outside r25- We have already seen that 
this de cline is still shallo w compared to the decline in 
EsFR. ILerov et al.l (|2008l) found that the CO and SFR 
scale lengths were comparable inside r25, so we expect 
that the decline in H I is also shallower than the decline 
in CO in the inner disk. This suggests that the total gas 
(H I -t- H2) distribution may follow a broken exponential 
rather than a single profile. 

In Figure [13] we compare individual line-of-sight mea- 
surements between the optical and outer disks. Plot pa- 
rameters largely match Figure [3 though now the outer 
disk data appear as unfilled contours (the median and 
estimated scatter are still black points with error bars). 
Filled contours sho w data from within the optical disk 
(jBigielet al.l 120081 their Figures 8 and 12 with shght 
changes to methodology: we assign equal weight to each 
galaxy, use a larger bin width, and define the contours 
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Fig. 11. — FUV-to-Hl ratio (SEE, y-axis) as a function of (logio of) Toonire's Q parameter (x-axis) in the outer disks of spiral (left) 
and dwarf (right) galaxies. Details of the plot are as in the previous figures. Q is calculated under the assumption of a fiat rotation curve 
and (Thi = 11 km s~^ ITamburro et ani2009l) . Almost all of the data are formally stable (Q > 1, right of the dashed line). There is no 
clear trend in the EUV-to-H I ratio (SEE) as a function of Q at ~ 750 pc scales, as one might expect if Q was the main criterion for a 
step-function type star formation threshold. 
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ElG. 12. — Scale lengths of EgpR (EUV + IR) measured within the optical disk by ILerov et al.l II2008I ) as a function of our outer disk 
SspR (EUV) scalelengths for the overlap in samples. The dashed line indicates equality and the dotted lines show ±0.1 r25. The two 
quantities agree reasonably well, with only the dwarf galaxies IC 2574 and Ho I showing a sharp truncation in Sjji as well as in Sgpjj . 



via the fraction of data enclosed rather than absolute 
counts). Tables containing the distribution and sampling 
data used to construct Figure [13] (optical and outer disk 
data for the spiral and the dwarf sample) are available 
as electronic supplemental material. 

The outer disks in Figure [T2] largely extend the distri- 
bution found within the optical radius to lower SFR and 
gas surface densities. There may be a small discontinuity 
between the two distributions along the y-axis due to the 



inclusion of an IR-based extinction correction inside the 
optical radius. What is striking is that the outer disk 
data lie overwhelmingly in a different part of parameter 
space from data inside the optical disk. In both samples, 
we see a smooth trend extend to very long depletion times 
and low gas columns that are almost totally absent inside 
the optical disks. Also in contrast to the optical disks, 
a clear trend relating H I and SFR emerges in both sub- 
samples. Inside r25, particularly within 0.5 x r25, the re- 
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Fig. 13. — Star formation and gas from inner to outer disks (left: spirals, right: dwarf galaxies). Filled contours show the pixel-by-pixel 
distribution of SgpR as a function of Sgas inside r25 (filled contours) for a subset of our sample (Bigiel et al. 2008). Empty contours show 
outer disk data (see Figure [Sjl with black filled circles indicating median and scatter. Plot parameters are as in Figure |8| In spirals, the 
steep relation between SsFR ^^nd Egas seen inside r25 becomes more shallow at large radii, leading to relatively modest variations in the 
median SFE (depletion time) across the outer disks (though the scatter is large). For both spiral and dwarf galaxies, data from inside r25 
appear to more or less continuously extend into the outer disk. Note that the slight offset in SsFR between the outer and inner disk data 
is due to different methodology in deriving SFRs in the two regimes (see text for details). These data are also available as supplemental 
online material. 



lationship between H I and Ssfr is weak or nonexistent; 
inste ad H2 is clearly correlated with Ssfr (Bigiel et al. 
[2OO8I) . 

There appears to be a 'forbidden region' in Figure [13] 
at relatively high Shi (~ 3-10 Mq pc~^ or log(SHi) be- 
tween ~0.5 and 1) and low Ssfr- This lack of high- 
H I low-FUV gas combines with the upper envelope of 
■''Dep ~ 10^*^ yr in the outer disks and the turn towards a 
fixed TDep in the H2 dominated (inner) parts of galaxies 
to create a combined distribution with an 'S-like' shape 
(though note that there is some gas at low columns with 
T'Dep ^ 10^° yr inside the optical disks of spirals). We 
will interpret this combined distribution in Section |4l 

The dwarf distributions lack the upper turn to the right 
in the 'S-shape:' we do not observe a significant amount 
of data at high Sgas, at least partially because we lack in- 
formation on the amount of H2 in these systems. Instead 
we observe a soft upper limit to the H I surface density of 
~ 10 Mq pc~^ at our resolution of 15" (the 'saturation' 
discussed by [Bigiel et al.l |2008[ ) . Figure [T3] shows that 
this pile-up at the saturation extends into the outer disks 
of dwarf galaxies, though the SFR for Ehi ~ 10 Mq pc~^ 
outside r25 is lower than for the same gas inside r25, on 
average. Unless H2 represents the majority of the ISM by 
several times in dwarf galaxies, it seems safe to conclude 
that as in spirals star formation is more efficient inside 
the optical disks of dwarf galaxies than at large radii. 

Finally, in Figure [141 we look at how the SFE varies 
as a function of galactocentric radius from to 2 x r25 ■ 
We plot our data at > r25 (right of the vertical dashed 



line) and measurements bv lLerov et al.l (|2008D inside r25 
(their Figure 1); we also indicate their fits to the SFE 
as a function of radius for each sample. As in the case 
of the SFR vs. gas plot (Figure [TH) , some discontinuity 
reflecting the change in methodology is seen around r25 
(we illustrate the effects of neglecting internal extinction 
within r25 by the black set of contours, which show an 
extension of the red filled contour of the outer disk data 
towards smaller radii). Our approach here is not appro- 
priate to investigating the detailed behavior around r25, 
nor is this our goal; we wish to broadly compare the inner 
and outer disk regimes. 

Again, we see a distinct contrast between the inner 
and outer disk. The SFEs overall are lower, the scatter 
in SFE appears larger, and although the first few points 
in both samples are approximately con sistent with the 
decline measured by iLerov et"aLl (|2008D , the overall de- 
cline in SFE with radius across the range 1-2 x 7^25 is 
markedly more shallow than between and 1 x r25. 

4. DISCUSSION 

4.1. Broad Structure of Star Formation In Outer Disks 

Our most basic conclusion is that star formation in 
outer disks is extremely inefficient compared to star for- 
mation inside the optical disks. H I depletion times for 
most of our data are ~ 10^^ yr and it is rare to find re- 
gions with local depletion times < 10^" yr. These values 
agree with previous measurements near r 25 and obser- 
vations of low surface brightness galaxies (jWvder et al.l 
,2009.) . Thus present day star formation requires many 
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Fig. 14. — SFE (i.e., SFR/gas or FUV/Hl) as a function of galactocentric radius from to 2 X r25 in spiral (left)and dwarf (right) 
galaxies. The vertica l dashed line indicates r25. Contours and black filled circles right of this line are as in Figure |9] Data inside r2s 
are from ILerov et al.l pOOS) and are derived accounting for internal extinction and the presence of molecular gas (meaning that some 
discontinuity at r25 may be expected). We illustrate the effects of neglecting these two factors by the black set of contours, w hich show an 
extension of the outer disk data (the red contour) towards smaller radii. The dashed lines show fits to the SFE inside r25 from ILerov et ahl 
112003 ). The outer disk data are roughly consistent with these trends inside ~ 1.5 r25 but at large radii the decline in the SFE is more 
shallow than it is inside the optical disk. 



Hubble times (~ 10^" yr) to devour the existing gas reser- 
voir. Another way to look at this is the low integrated 
SFR of outer galaxy disks: on average, outer disks con- 
tribute only about 10% to the total SFR of a galaxy. 
The relative lack of importance of in situ star formation 
means that the massive extended gas distributions ob- 
served in many nearby galaxies can be long-lived and may 
be viewed as a potential source of fuel for inner disk star 
formation. The short depletion times in inner disks {^ 
2x10^ vr. lBigiel et al.l2 "008':'Lerov et al."2008') imply that 
such a source is required ( Shlosman, Frank & BegelmanI 
[1989; Blitd 119961 iBauermeister et ah 2009), though the 
presence an d importance of radial gas flows are still de- 
bated (e.g ., IVollmer fc BeckertI [200a IWong et al.l [20p: 
[Peea[2009h . 

Our second simple conclusion is that both H I column 
and FUV intensity decline systematically at large radii. 
For H I, this represents a contrast with the inner disk, 
where the azimuthally averaged H I surface density tends 
to be flat or even increase with increasing radius. As a 
result, outer disks are home to column densities seldom 
found inside the optical radius; in the regime we study 
in spirals, column densities < 3 Mq pc~^ account for the 
majority of data. For FUV, the decline that we observe 
outside r25 appears approximately consistent with the 
decline in SFR inferred from com bined IR and UV pro- 
files inside ^25 (jLerov et al.ll2008D . The rate of decline in 
H I and FUV intensity for a given galaxy appears to be 
correlated but not identical; parameterizing the decline 
with exponential scale lengths, we find the FUV inten- 
sity to decline with a scale length approximately half 
that of the H I. With the SFR declining faster than the 
gas reservoir, star formation becomes less efficient with 
increasing radius, though this decline is actually milder 
than that observed inside r25. 

Despite the mismatch in scale lengths, we find a fairly 



good overall relationship between H I and FUV (mea- 
sured via direct scatter plot or rank correlation). This is 
surprising given the extremely poor correlation between 
H I and SFR (trace d mainly by IR emission) found inside 
star-forming dis ks (Wong & Blit^'2002^ ' Kennicutt et al.l 
I2007t ILerov et al. 2008;"Bigicl ct al. 200^) and gives a 
strong hint that different physics governs the formation 
of star-forming clouds at lar ge radii. This agrees qual- 
itatively with the results of iHunter et al.l ()2010D . who 
found that the integrated Hl-richness of dwarf galaxies 
appears related to outer disk star formation in the sense 
that Hl-poor dwarfs preferentially show suppressed star 
formation. 

4.2. Star Formation Thresholds 

Although we do find the SFE to depend on both gas 
column and radius, we do not find clear evidence for a 
threshold in either quantity. This agrees with numerous 
recent studies that have failed to detect s uch a feature 
using UV data (e.g., iBpissier et al. 2007^ jThilker et alj 
l2005tlllmrt"er et al.ll2010l : IGil de Paz eral1l2005LI2007bD . 



but is in c ontrast to the rad i al cutoffs found in Ha emis - 
sion (e.g., iKennicutti [19891: iMartin fc Kennicutti [200l . 
The main new contribution of our work is to look for the 
signature of such a feature by plotting the SFE of indi- 
vidual lines of sight as a function of H I column density. 
Doing so, we observe a continuous dependence of SFE 
on H I column with large scatter rather than an abrupt 
step at a given H I column density. 

What does the lack of a step function in FUV/H I as 
a function of H I mean? Our angular resolution (15") 
translates to ^^ 750 pc for a typical spiral in our sample 
(d ~ 10 Mpc). A reductionist view of our observations 
would imply that if there is a critical column density that 
must be achieved on these scales in order for star forma- 
tion to occur, then this column density is not maintained 
over the same timescales and area as the FUV emission 
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from the young stars that resuh. Our best estimate is 
that this is not due to internal extinction or unaccounted 
H2, though in both cases this represents an assumption. 
It may be a timescale effect, with FUV emission surviving 
its parent gas. FUV emission at low but significant levels 
may be found from stellar populations as old as 100 Myr, 
during which a population with a 10 km s~^ dispersion (a 
plausible dispersion for young stars and roughly the ob- 
served dispersion in the gas) might move ~ 1 kpc. This 
would allow FUV emission from one resolution element 
in our line of sight analysis to diffuse to nearby regions 
and the gas distribution to evolve significantly. 

Another explanation is that star formation has only a 
stochastic relation to kiloparsec scale conditions. In the 
inner parts of galaxies, star forming clouds are ^ 50 pc 
in size, a factor of 10-20 times smaller than one of our 
resolution elements. A similar contrast holds in terms of 
gas mass: even at a relatively modest Ehi ~ 1 Mq pc~^, 
one of our resolution elements contains '^ 5 x 10^ Mq 
of H I. This is enough material to make roughly 10 low- 
mass star-forming regions the mass of the Taurus molec- 
ular complex — perhaps a reasonable analog to outer 
disk star forming regions given the low local SFRs and 
sparsity of high mass stars. Detailed GALEX studies 
of outer disks support this picture, finding star forma- 
tion to be patchy, composed of locally confined 'FUV 
knots' (|Gil de Paz et al.ll2007at iThilker et ai][2007l ) that 
may match theoretical exp ectations for this regime (e.g., 
lEhnegreen fc H unter 2006\ 

An even simpler explanation can be obtained by real- 
izing that local conditions clearly influence the amount 
and efficiency of star formation in systematic, measur- 
able ways. However, the convolution of a turbulent ISM, 
random sampling of star-forming clouds along their evo- 
lutionary sequence, and the sensitivity of star formation 
tracers to a range of stellar ages means that star for- 
mation 'thresholds' will always be obse rved as contin - 
uous trends. The physi c s disc ussed bv iSchavd (I2004D 
lElmegreen fc Parravand (|1994[ ). and iKennicuttI (|19890 
may drive the trends that we observe, but we expect it 
is unlikely that any proposed threshold will yield a clear 
step function on kpc scales (and measurements on much 
smaller scales risk returning the trivial result that stars 
form in dense, bound clouds rather than assess where 
stars can and cannot form on galactic scales). 

4.3. What Regulates Star Formation In Outer Disks? 

We have seen outer disks to be distinct from the in- 
ner parts of spirals in several ways: SFEs are low, the 
gradient in SEE with radius is comparatively weak, and 
there is a clear spatial correlation between H I and SEE. 
We have also found the outer disks of dwarf and spiral 
galaxies to be similar in many ways. All of these, we 
argue, point to H I column as the regulating quantity 
for star formation in outer disks. There are two natural 
reasons for this: first, H I represents the ultimate fuel 
for star formation and unlike in inner galaxy disks the 
availability of H I varies dramatically across outer galaxy 
disks. Second, inasmuch as the outer parts of galaxies 
are actually organized into disks, the H I column will 
usually represent the dominant (baryonic) mass compo- 
nent, meaning that the volume density of gas (which de- 
pends critically on the stellar surface density inside r25) 
depends mostly on Ehi- Typical rank correlation coeffi- 



cients across the outer disks of our galaxy sample show 
a comparable degree of local correlation (r « 0.4) be- 
tween the SEE and H I column and between SEE and 
radius. This can be (at least partly) understood as both 
quantities are strongly anti-correlated (i.e., are not inde- 
pendent; compare Figure [2]). The correlation coefficient 
between SEE and Toomre-Q (r < 0.1) implies little or 
no local correlation between these quantities. 

To first order a line of sight with a given H I column 
appears to form stars in a way that is largely independent 
of whether it lies in the outer disk of a dwarf or spiral 
galaxy. There are second order differences that point 
to the importance of other factors, however. The most 
straightforward of these is the large scatter in Ssfr at a 
given H I column in spirals. Our estimate of the intrin- 
sic scatter in Esfr near Ej^/ ^ 1 M© pc^^ is w 1 dex. 
Some of this might be expected from the low values in- 
volved: the median SER for the lowest bins corresponds 
to forming only a few hundred solar masses every 100 
Myr. However, over the same range dwarf galaxies show 
significantly less scatter; with such a small sample of 
dwarfs it is hard to know whether the difference is sig- 
nificant, but if it is then an easy interpretation is that 
the larger scatter in spirals reflects the wider range of 
environments found there (e.g., our sample includes spi- 
rals with strong metallicity gradients like M 101 and such 
with shallow gradients like M 51 and M 83). 

What can we say about the efficiency of the regulation 
of star formation in outer disks? In the inner parts of 
galaxies stars form out of molecular clouds with a deple- 
tion time of a few times 10^ yr. Across both dwarfs and 
spirals we find H I depletion times in outer galaxies two 
orders of magnitude higher, '^ 10^^ yr. If an outer disk 
star-forming cloud resembles its inner disk counterpart 
(a big assumption), then only ~ 1% of the gas in outer 
galaxy disks is actually in star forming clouds (or about 
~ 5,000 M0 in one of our typical resolution elements). 
Similarly, if we take the orbital time to be ~ 0.5 Gyr 
(roughly appropriate for the Milky Way around r25) then 
only ^ 0.1% of the gas is converted into stars per or- 
bital time, two o rders of m agnitude below the galaxy- 
averaged value of IKennicuttI (1998 ). Assuming a typical 
scale height of ~ 500 pc, average volume densities in the 
outer disks will be ~ 0.1 cm~^, implying a disk free- 
faU time of ~ 10^ vr. iKruinholz fc McKed pOOl and 
iKrumholz fc Thompson! (j2007f ) argued that star-forming 
structures on many scales convert '^ 1% of their gas to 
stars per free-fall time. Outer galaxy disks apparently 
have about a tenth that efficiency on average. 

4.4. Towards a Complete Star Formation 'Law' 

Finally, we comment on the apparently 'baroque' na- 
ture of the 'star formation law' plot seen in Figure 1131 
The plot shows a complex star formation law, with 
at least three (somewhat) distinct regimes observed — 
outer disks, the H I-dominated parts of inner disks, and 
the H2-dominatcd parts of inner disks. Another regime 
is well-established: starburst galaxies are known to have 
a steeper power law index than the H2- dominated parts 
of spirals in Figure [131 so that the plot should 'turn up' 
again at Sgas ^ 100 Mq pc~^ (compare discussion in 
iBigielet alTl2008D . Different parts of the 'S-shape' (a 
'W-shape' with the starbursts) represent physically truly 
distinct regimes. We argue that this kind of behavior — 
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a series of regime- dependent relations that combine into 
a complex distribution — is actually exactly what one 
would expect if the ability to form stars or star-forming 
clouds has a significant environmental dependence and 
multiple environments are combined on a single plot. 

5. SUMMARY 

Using high resolution ('^ 15"), sensitive, wide field- 
of-view data we have studied the relationship between 
H I gas and star formation traced by UV light in the 
outer disks (1-2 x r25) of 17 spiral and 5 dwarf galaxies. 
We expect that H I dominates the ISM in this regime 
and that FUV reflects the distribution of recently formed 
stars without large biases due to internal extinction. 

We find that: 

1. Despite widespread star formation the gas deple- 
tion time is very low in outer galaxy disks, almost 
always longer than 10"'^'^ yr and typically ~ 10^^ yr. 
Star formation at the present rate will thus not ap- 
preciably deplete the gas supply in the outer disks 
over timescales shorter than many Hubble times. 

2. FUV emission in outer disks declines approxi- 
mately exponentially with a scale length compa- 
rable to that found for star formation in the inner 
parts of galaxies. H I gas also usually shows a sys- 
tematic decline, but a shallower one than the FUV. 
The H I scale length that we derive for outer disks 
is typically twice the FUV scale length. As a result, 
the ratio of H I to FUV and thus the gas (H I) de- 
pletion time, increases with radius. However, this 
increase is weaker than that observed inside r25 , so 
that overall the depletion time as a function of ra- 
dius becomes 'flatter' as one moves out in a galaxy. 

3. Unlike in the inner parts of galaxies, there is a fairly 
clear relationship between star formation and H I 
in outer galaxy disks. This is clearest near the op- 
tical radius, where the H I and FUV are strongly 
correlated (rank correlation '^ 0.7). Extinction 
and the presence of H2 wash out this correlation 
at radii < r25 and the strength of the correlation 
also decreases with increasing radius (r > r25), 
perhaps due to the increasing scatter/stochasticity 
in star formation at these radii. A fit of deple- 
tion time vs. Shi to binned outer disk data yields 
TDep oc (Shi)°'^, or SsFR oc (Ehi)^'^ for both spiral 
and dwarf galaxies. 



6. The relationship between Ssfr and Sgas for in- 
dividual lin e s of s ight extends that measured by 
iBigiel et"al] ()2008D into a new regime of low Esfr 
and low Sgas more or less continuously. In this 



regime, a line of fixed td, 



cp 



10 yr forms an 
approximate upper envelope to the distribution. 
Lines of sight with high H I columns (Sgas ~ 3- 
10 Mq pc~^) but little or no FUV emission are 
rare. These conclusions apply to both dwarf and 
spiral galaxies. 
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At a given H I column density, one finds similar 
FUV intensities in the outer disks of dwarf and 
spiral galaxies. The main difference between the 
two samples is the distribution of H I column den- 
sities at large radius. Dwarf galaxies tend to have 
(relatively) more high-column H I than spirals. In 
both samples, low columns become more common 
with increasing radius. 

Examining the relationship for individual lines of 
sight, we do not find clear evidence for a step- 
function type star formation threshold in either 
radius or H I column, though the SEE is a weak 
function of both quantities. 
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